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Abstract 

We study the stability of static black holes in the third order Lovelock theory. We derive a 
master equation for tensor perturbations. Using the master equation, we analyze the stability of 
Lovelock black holes mainly in seven and eight dimensions. We find there are cases where the 
linear analysis breaks down. If we restrict ourselves to the regime where the linear analysis is 
legitimate, black holes are always stable in seven dimensions. However, in eight dimensions, there 
exists a critical mass below which black holes are unstable. Combining our result in the third order 
Lovelock theory with the previous one in Einstein-Gauss-Bonnet theory, we conjecture that small 
black holes are unstable in any dimensions. The instability found in this paper will be important 
for the analysis of black holes at the LHC. 
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I. INTRODUCTION 



It is believed that string theory is a promising candidate for the theory of everything. 
Remarkably, string theory can be formulated only in ten dimensions. Apparently, it is 
necessary to reconcile this prediction with our real world by compactifying extra-dimensions 
or by considering braneworld. Recently, the idea of large extra-dimensions in the context 
of the braneworld has been advocated IlL Intriguingly, in the presence of the large 

extra-dimensions, black holes could be created at the TeV scale Qj. Hence, the stability of 
higher dimensional black holes becomes important since these black holes could be produced 
at the Large Hadron Collider (LHC) if the spacetime has larger than six dimensions. 

The stability of higher dimensional black holes has been an active topic since the seminal 
papers by Kodama and Ishibashi Q]. There are at least two directions to be pursued. One 
is to study various black holes in Einstein theory. This direction is necessary because black 
holes produced at the LHC are expected to be charged or rotating. A numerical study of 
charged black holes has been done To investigate the stability of rotating black holes, 
a group theoretical method is developed pj. The method is used to study the stability of 
squashed black holes js, [^j and 5-dimensional rotating black holes 
rotating black holes in more than 5-dimensions is also studied 
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10] . The stability of 
13, UJ1. The other 



direction is to consider the stability of black holes in more general gravitational theories. 
This direction is also important because black holes are produced at the Planck scale where 
Einstein theory would be no longer valid. In fact, it is known that Einstein theory is merely a 
low energy limit of string theory [lsj] . In string theory, there are higher curvature corrections 



in addition to Einstein-Hilbert term [15[. Thus, it is natural to extend gravitational theory 
into those with higher power of curvature in higher dimensions. It is Lovelock theory that 
belongs to such class of theories [if], In Lovelock theory, it is known that there exist 



static black hole solutions [18J. Hence, it is natural to suppose black holes produced at the 
LHC are of this type [lj]]. Thus, it is important to study the stability of these Lovelock 
black holes. 

In the case of the second order Lovelock theory, the so-called Einstein-Gauss-Bonnet the- 
ory, the stability analysis under tensor perturbations has been performed 20| (see also an 



earlier work 



2 ll|). The analysis has been also extended to the scalar and vector perturba- 



tions by the same group 22]]. They have shown there exists the scalar mode instability in 
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five dimensions, the tensor mode instability in six dimensions, and no instability in other 
dimensions. However, in the dimensions higher than six, we need to incorporate higher 
order Lovelock terms. Indeed, when we consider black holes at the LHC, it is important to 



consider more higher order Lovelock terms 



231 ] . Hence, in this paper, we study the stability 



of black holes in the third order Lovelock theory. In this context, our discussion of the 
stability of black holes in seven and eight dimensions becomes complete. In this paper, we 
will restrict ourselves to the tensor mode analysis as a first step. Even in this case, we find 
the instability of small black holes in eight dimensions. This suggests the instability of small 
black holes is generic in Lovelock theory. 

The organization of the paper is as follows. In Section [III we review Lovelock theory and 
introduce Lovelock black hole solutions. In Section IIII[ we present a master equation for 
tensor perturbations in the background of Lovelock black hole. In Section IIV| we clarify 
the conditions of applicability of the linear analysis and the stability. Then, we discuss the 
stability of Lovelock black holes for various cases. In particular, we prove the stability of 
black holes under tensor perturbation in seven dimensions and find the instability of small 
black holes in eight dimensions. Finally, we summarize our results in Section |Vj 

II. LOVELOCK BLACK HOLES 

In this section, we present the third order Lovelock theory. In the spacetime dimensions 
D < 6, the third order Lovelock term is not relevant. While, for D > 9, we need to 
incorporate the fourth order Lovelock term, which is beyond the scope of this paper. Hence, 
our main objective in this paper is the third order Lovelock theory in seven and eight 
dimensions, although we mention black holes in other dimensions. 



In 16], the most general symmetric, divergence free rank (1,1) tensor is constructed out 
of a metric and its first and second derivatives. The corresponding Lagrangian is given by 

l(D-l)/2] 

L = ^ ] c n C m , , (1) 

m=0 

where C m is defined by 

r a m bm p c\d\ p 

J -"m — 2m°cid 1 -c m d m It a, 1 bi - rt a m 6 m • 

Here, [z] means the largest integer satisfying an inequality [z] <z, c n is an arbitrary constant 
and fi^di-'-^dm * s ^ ne generalized totally antisymmetric Kronecker delta. 



In this paper, we consider up to the third order term in the Lagrangian ([I]). From now 
on, we set Co = — 2A, c\ — 1, c 2 — a/2, c 3 = f3/3, and c m = (for 4 < m). Then, we can get 
the following equation from this Lagrangian: 

= Q b = A5 a b + G^ b + aG^ b + f3G^ b , (3) 

where G^ b = R a b — \R5 a b is the Einstein tensor. The second order Lovelock tensor, i.e., 
the so-called Gauss-Bonnet tensor, G^ b is given by 

/^f(2)6 p de p cb r> p c p db OE>cpb,ppb A fc /* f/H 

_ -"-co -Kde ~ ^tid tica ~ ^ti a tic + KK a — —0 a L, 2 ■ (4 J 

The third order Lovelock tensor Ga 3 ^ b reads 

/^»(3)6 p2 p b /(pfeprfpci pbp e / E> ccf ^ DD bd-Q c 

(jr a ' — K K a — 4/t a /t c ild + /t a ilcrf /t e / — 4:KK ac Kd 

I op ba! p c/pe L op fedpcpe zip bd p cop e/ 
+ o-n-ac -Krfe -"-/ + °K ac K e Kd — *iK ac K e f ridg 

/IDpcpfe I o p b p ce r> d . o p c p d p b ,1 p ip c / p ^ e 

— 4itrt a /t c + 8/t c rt a rf /t e + Q/t a it c Kd — QKc Kd e K a f 

i o p p ffe p 6c /) p de p b/ p c , a p c p b/ p de , r> p e/ p be p g/i 
+ lKK ac Kde — ^ti ac Kde tij + 4r£ a rf^e it/c + ZK ac K g h K e f 

opaiepfecp/iopcpfcepa! Q D p bg p /c Jf & /* /"cN 

— o-rt ac -rtrf/ K e + oK a K c d K e — QK ac Kdf K ge — -0 a L,i . [0) 

In the above tensors, the Lagrangians 

C 2 = R 2 - 4R c d Rd c + Red ef Re f cd (6) 

and 

£3 = R 3 — l2RR c d Rd c + 16R c d Rd e R e c + 2R C d e R e f 9h R g h cd 

+8R cd ef Re g ch R hf d9 + 3RR cd ef Re f cd - 2AR cd ef R eg cd R f 9 + 24R cd ef R e c R f d . (7) 

have appeared. 

As is shown in [la ], there exist static black hole solutions. The line element can be 
expressed as 

ds 2 = -f{r)dt 2 + y^y^r 2 + r%jdx l dx j , (8) 

where 7^ is the metric of n = D — 2-dimensional manifold with constant curvature k = 1, 
or — 1. The function f(r) is determined by Eq.(j3]). It is convenient to define a new variable 

ip as 

f(r) = K-r 2 ^(r) . (9) 
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FIG. 1: The intersection between the solid curve and thin horizontal line determines the solution 
ip = ij){r) for the case n = 4. Apparently, the infinity r = oo corresponds to ?/> = 0. The intersection 
between solid and dashed curve gives a horizon r#. 



Then, the solution can be found from the following implicit equation 



f (n - l)(n - 2)(n - 3)(n - 4)^ 3 + £(n _!)(„_ 2 )^ 2 + rj, - = -£L 

3 2 n(n + 1) r n+1 



, (10) 



where \i is a constant of integration related to ADM mass 24j : 

o u7r (n+l)/2 

M = ^- — -, (11) 

r((n+l)/2) 1 ] 

where we used a unit 167rG = 1. 

In this paper, we will concentrate on the asymptotically flat spherically symmetric solu- 
tions, i.e., A = and k = 1. We consider only positive mass black holes \x > 0. In addition 
to these, we impose extra assumptions a > and j3 > 0, for simplicity. 

In the case of n = 3, the theory is reduced to Einstein-Gauss-Bonnet theory. In this case, 
we can solve Eq. ffTOl) explicitly 



-1 ± 



4a fi 



2a 



(12) 



The upper branch leads to the asymptotically flat solution 

.2 



f(r) 



r 

l + 2~a 



1 + 



4a fi 



(13) 
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FIG. 2: The same method as n = 4 case is illustrated for n = 5 case. 



In the case of n = 4, we have a similar expression. In Fig{TJ a graphical method is also 
explained for this case. In other dimensions, analytic formula would not be useful. Hence, 
we illustrated a graphical method for n = 5 in FigJ2 We should note that in order to 
have asymptotically flat solutions we have to restrict solution of Eq. (fT0l) to the positive one. 
Then, as one can see easily from Eq. (fT0|) . ifj goes as l/r n+1 in the asymptotic region. Hence, 
f(r) gives the asymptotically flat metric. 

The horizon radius of the asymptotically flat solution is characterized by /(V#) = 0. 
From ([9]), we have a relation ipn = ^{fn) = l/ r /f- Using this relation and ( JTU1) . we obtain 
an algebraic equation 

|(n - l)(n - 2)(n - 3)(n - 4)^ + |(n - l)(n - 2ty| + ^ = rf +1)/2 . (14) 

This determines and hence r^. In seven dimensions, one can calculate ipn analytically. 
In fact, substituting n = 5 into (1T41) . we obtain a simple equation 

(8/3 - /i)^fr + 6«Vh + 1 = ■ (15) 
Then, if ji > 8/3 — 9a 2 , the solution is given by 



3a + v/9a 2 - 8/3 + fi 
= ^ • (16) 

In other dimensions, we do not have analytic formulas. However, it is easy to obtain solutions 
of (fl4j) numerically. The method to solve Eq. (|T4|) is illustrated in FigJUand Figj2]in the case 
of n = 4 and n — 5, respectively. We note that ip moves from to ipu- We also notice that 
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the larger ipn gives the smaller \x in this case. This is clear from Figfj] and Figj2j This is 
important for later discussion of the stability. 

From the metric (|S]), we can calculate the Kretschmann scalar R a b a iR abcd as 

R abcd R abcd = f" 2 + + 2n{n - , (17) 

where the prime denotes a derivative with respect to the coordinate r. Thus, the solution has 
curvature singularities at r = or other point where derivative of / diverges. The behavior 
of a solution can be understood from Eq. ljTOl) . In fact, the solution satisfy the conditions 
|1 — f \ < oo, \f'\ < oo and \f"\ < oo except for r = 0. Therefore, a curvature singularity 
exists only at r = 0. 

Thus, there is an asymptotically flat solution with a horizon hiding a singularity at r = 0. 
Hence, this solution describes a black hole with the mass M defined in ffTTl) . 



III. MASTER EQUATION FOR TENSOR PERTURBATIONS 



Understanding of perturbed black holes is important for the analysis of black holes at the 
LHC. Of course, the stability of black holes is a prerequisite for this. In the case of higher 
dimensional black holes, the stability is not guaranteed in contrast to the 4-dimensional 
cases. To prove the stability, it is convenient to decompose the metric under the symmetry 
of n-dimensional symmetric space. There are scalar, vector, and tensor modes. In this 
paper, we will concentrate on the tensor perturbations as a first step. 

We consider tensor perturbations around the solution ([H]) 

Sg ab = , Sg ai = , Sg^ = r 2 (f}(t,r)h i:j (x l ) , (18) 

where a,b = (t, r) and 0(t, r) represents the dynamical degrees of freedom. Here, hij are 
defined by 



V k V k h 



k't'ij 



IK 



13 i 



0, 7% 



0. 



(19) 



Here, V* denotes a covariant derivative with respect to 7^. Here, the eigenvalue is given by 
7 = — £(£ + n — 1] + 2, (& — 2, 3, 4 • • • ) for k = 1 and negative real number for n = —1, 0. 

As is shown in 20|, tensor perturbations around the solution ([8]) in Einstein- Gauss-Bonnet 
theory can be calculated as 



SG) 



(i)i 



(20) 



and 



5G) 



(2)3 



/¥'; 

n-2 

2r 3 

7 - 2k 
2r 4 



n-2 

w 



{-r/ + (n-3)(«-/)} 



{(„_3)[( n _2)(/-«)/- r «/] 

+r 2 (f' 2 + /"/) + (3n - 7)r/'/} 
[ r y + 2 ( n _ 3 )r/ + (ra - 3)(n - 4)(/ - «)] 



(21) 



After a long calculation, we also obtain the third order contribution 



5G) 



(3)i 



t»\ 1 



, 1 



-0_(n-2)(n-3)(n-4) 



(/ - K){2r 2 //" + (5n - 21)r// - (n - 5)«r/ 



+ (n - 5)(n - 4)/(/ - k) + 4r 2 /' 2 \ + 



,2 r'2 



+^^_7(n-3)(n-4) 



(/ - «){(n - 5)(n - 6)(/ - k) + 4(n - 5)r/ + 2r 2 /"} + 2r 2 /' 2 
Given the formula (1201) . (12 ip and (1221) . the equation can be written as 



5G 



aSG) 



(2)3 



PSG\ 



(3)3 



. 



hi .(22) 



(23) 



Separating the variable 0(r, t) = x{ r ) eWt i we can deduce the equations to the following form: 

(24) 



(n — 2)r /i 



where 



/i(r) = r n " 2 + a(n - 2)r n ~ 4 {-r/' + (n - 3)(k - /)} 

+ /3{n - 2){n - 3){n - 4)r n " 6 (fi: - /) {-2r/ + (n - 5)(« - /)} . (25) 

The sign of h and its derivative dh/dr is crucial for the stability analysis. 

For later purpose, it is useful to express the function h not by the coordinate r but by ip 
defined in ( TTUj) . From Eq. ({T0l) . we can deduce a relation 



[P(n - 4)(n - 3)(n - 2)(n - l)ip 2 + a(n - 2)(n - 1)V> + l] ^' 

~(n - l)(n - 2)(n - 3)(n - 4)^ 3 + -(n - l)(n - 2)^ 2 + $ 



(26) 



Using Eq. (fTO|) and the above relation (I26I) . we can express (I25I) as a function of ip: 



- r «-2 hktl 

6[/3(ri-4)(n-3)(n-2)(™-l)^ 2 + a(n-2)(n-l)^ + l]' v ; 



where 



L(x) = 2(n-5)(fi-4f(fi-3) 2 (n-2) 2 (n-l)/3V 
+4(n - 5)(n - 4)(n - 3)(n - 2) 2 (n - l)a/3x 3 
+3(n - 3)(n - 2) {(n - 2)(n - l)a 2 - 8(n - 4)/3} 
+6(n-3)(n-2)ax + 6 . (28) 



x 2 



Similarly, dh/dr is given by 

^ = r n-3 (29) 

rfr 36 \j3(n - 4)(n - 3)(n - 2)(n - l)^ 2 + a(n - 2)(n - 1)^ + l] 3 ' 

where 

K(x) = 4(n - 8)(n - 5)(n - 4) 4 (n - 3) 4 (n - 2) 4 (n - l) 3 /5 4 x 8 

+ 16(n - 8)(n - 5)(n - 4) 3 (n - 3) 3 (n - 2) 4 (n - l) 3 a/?V 

+2(n - 4) 2 (n - 3) 2 (n - 2)\n - lf{{n - 2){n - l)(16n 2 - 163n + 469)a 2 

-4(n - 4)(n - 3)(n + 13) (2n - 7)(3}(3 2 x 6 
+6(n - 4)(n - 3)(n - 2) 3 (n - l) 2 {(n - 2)(n - l)(5n 2 - 40n + 99)a 2 

+2(n - 4)(n - 3)(n 2 - 43n + i72)/3}a/3a; 5 
+3(n - 3)(n - 2) 2 (n - l){3(n - 5)(n - 2) 2 (n - 1) V 

+2(n - 4)(n - 2)(n - l)(13n 2 - 104n + 315)a 2 /3 
+8(n - 4) 2 (n - 3)(n 2 - 8n + 45)/3 2 }x 4 
+12(n - 3)(n - 2) 2 (n - l){3(n - 5)(n - 2)(n - l)a 2 

+2(n - 4)(7n 2 - 21n + 80)/5}ax 3 
+ 18(n - 2){(n - 2)(n - l)(5ra 2 - 25n + 42)ct 2 

+4(n - 4)(n - 3)(2n 2 - n + 9)p}x 2 
+ 108(n - 2)(n 2 - 3n + 4)ax + 36(n - 2) . (30) 

In the next section, we will study the stability of black holes using these formulas (1271) . 
([28]), ([29]), and (HP). 
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IV. STABILITY ANALYSIS 



Before doing the stability analysis, we will clarify the conditions for the stability of black 
holes in Lovelock theory. 

First of all, it should be noticed that we have to impose the condition 

h(r) > , (for r > r H ) ■ (31) 

This is necessary for the linear analysis to be applicable. If we choose f(r) to be asymptot- 
ically flat solution, one can easily confirm that /(r) ~ k — fi/r n+1 for large r. It is easy to 
see h(r) is positive in the asymptotic region. Therefore, the condition (I3~lj) means that the 
equation h(r) = has no solution in the region r > r#. In the case that there exists tq such 
that h(r ) = and r > rjj, we encounter a singularity. To see this, let us examine (|24|) 
around r . Using approximations h{r) ~ h'(r )(r — r ) = h'(r )y, f(r) = f(r ) and r = r , 
we can reduce Eq. (1241) into the following form: 

where c is a constant. The solution near tq is given by x ~ c i + c 2 logy, where c\ and C2 are 
constants of integration. The solution is singular at y = for generic perturbations. The 
above calculation means that if the condition (I3T|) is not satisfied, then the linear analysis 
break down at r^. The similar situation occurs even in cosmology with higher derivative 
terms [25]. 



When the condition (|3~T|) is fulfilled, we can introduce a new variable 

$( r ) = x{r)ryjh(r) . (33) 

Using \l/ and switching to the coordinate r*, defined by dr*/dr = 1/ f, we can rewrite Eq. (j24|) 

as 

^ + V(r{r*))V = -u 2 m = , (34) 



dr 



*2 



where 

V(r) 

is an effective potential. 



(2k — 7)/ din h 1 f d f f d ^\ 
(n - 2)r dr r \fh 



4r (4r^) < 35 ' 
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For discussing the stability, the " S- deformation" approach is very useful {5), Q. Let us 
define the operator 

As ~^ + V < 36 > 

acting on smooth functions defined on / = {r* H , 00). Then, ()3~4l is the eigenequation and E 
is eigenvalue of A. In this case, for any tp, we can find a smooth function S such that 

(tp,Atp) = J(\Dtp\ 2 + V\tp\ 2 )dr*, (37) 

where 

D = ^ + S, V = V + ff-S\ (38) 
dr* dr 



Following [20l | . we choose 5 as 

5 = -flL\n(rVh) . (39) 

Then, we obtain the formula 

f i* 00 1 2 J 1 i. 

4,) = / + (2k - 7) j( 5^-5-* . (40) 



Here, the point is that the second term of (1401 includes a factor 2ft — 7 > 0, but h does not 
include 7. Hence, by taking a sufficiently large 2k — 7, we can always make the second term 
dominant. 

Now, let us show the importance of the sign of d In h/dr. If d In h/dr > on I, the solution 
(EJ) is stable. This can be understood as follows. Note that 2k — 7 > 0, then we have V > 
for this case. That means (p,Ap) > for arbitrary p if din h/dr > on /. We choose, 
for example, tp as the lowest eigenstate, then we can conclude that the lowest eigenvalue E 
is positive. Thus, we proved the stability. The other way around, if din h/dr < at some 
point in J, the solution is unstable. To prove this, the inequality 

T^T " Eo (41) 

is useful. This inequality is correct for arbitrary p. If d In h/dr < at some point in 7, we 
can find </? such that 

\p\ 2 dlnh , . 

. — . — < . 42 

. [n-2)r dr v ; 
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In this case, ( 140]) is negative for sufficiently large 2k — 7. The inequality (141~]) implies Eq < 
and the solution has unstable modes. Thus, we can conclude that the solution is stable if 
and only if dlnh/dr > on I. 

To summarize, we need to check the sign of h and dh/dr outside the horizon r > Th to 
investigate the stability. Note that the condition h(r) > is necessary for the consistency of 
the analysis. To study black holes which do not satisfy this condition, we have to go beyond 
the linear analysis. 

Now, we are in a position to discuss the stability of Lovelock black holes. Here, we 
consider the most important cases, namely, asymptotically flat black holes with a > 0, 
j3 > 0, A = 0, /i > and k = 1. Notice that ip moves on < ip < ipn if the solution is 
asymptotically flat. Note that ipu is defined in (I14p . Of course, the stability analysis for 
other cases can be easily done by using our formula. 



A. a=0 and 0—0 case 

Let us start with black holes in Einstein theory, i.e., a = (3 = 0. For this case, the 
function h(r) defined by (125]) becomes h = r n ~ 2 . Apparently, h > and dh/dr > hold. 
Hence, black holes are stable. This is consistent with the known result j^]. 



B. a > and /3=0 case 



Next, we consider Einstein- Gauss-Bonnet theory 20[. From ( |27l) and ( !28l) . h(r) can be 
read off as 

= r , n _ 2 (n - 3)(w - 2)\n - l)aV + 2(n - 3)(n - 2)a^ + 2 

2[a(n-2)(n-l)^ + l] ' 

It is easy to see h is always positive because ip > and a > 0. Namely, the linear analysis 
is always applicable in Einstein-Gauss-Bonnet theory. 

Then, we can proceed to check the signature of dh/dr to study the stability of black 
holes. From (129!) and (1301) . we obtain 

^ = {n _ 2]r n-3 EM (44) 

dr v ; 4[(n-l)(n-2)a^ + l] 3 ' V ' 
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where 



F[<0] = (n- lf(n - 2f(n - 3)(n - 5)«V + 4(n - l) 2 (n - 2) 2 (n - 3)(n - 5)a 3 ^ 3 
+2(n - l)(n - 2)(5n 2 - 25ra + 42)aV + 12(n 2 - 3n + 4)a^ + 4 . 

Since both ijj and a are positive, dh/dr is always positive if n ^ 4. This means that black 
holes are stable under tensor purturbations other than in six dimensions. 
In six dimensions, dh/dr becomes 

dh = k(ajj) 
dr (6«V + 1) 3 ' 

where 

k( x ) = -54a; 4 - 36a; 3 + 66a; 2 + 24x + 1 . (46) 

Because the denominator of (1451) is always positive, the signature of dh/dr is determined 
by k(x). Note that k(x) is always positive in the region < x < |(— 1 + a/25 + 10a/6) 
and always negative in the region x > ~(— 1 + a/25 + IOa/6)- Since ^ lies in < ip < ipjj, 



dh/dr is positive when the inequality aipn < h(— 1 + v 25 + 10a/6) holds. Substituting the 
marginal value aipn = h{— 1 + a/25 + 10a/6) into Eq.( fI4l) . we get the marginal \i\ 



3a/6 (l + a/25 + lOv^) 

/imarginal = ~ ^ 3/2 ^ ~ 3-98247a 3/2 . (47) 



'25 



10a/6J 



It is obvious that black holes are stable if \i > 3.98247a 3 / 2 . 

Summarizing this case, black holes are always stable other than in six dimensions and 
those with fx larger than the critical value // marg inal = 3.98247a 3 / 2 are stable in six dimensions. 



This is consistent with the previous result 



20]. 



Thus, we have reproduced known results. This is an evidence of correctness of our master 
equation. Now, we can proceed to obtain new results. 



C. a=0 and (3 > case 

Let us consider a special case where the Gauss-Bonnet term accidentally vanishes. For 
this special case, we can give rather general results. 
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n=5 l(x) and k(x) 
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FIG. 3: The behavior of l(x) and k(x) in seven dimensions is shown. Both the lines cross the 
x-axis. Notice that l(x) crosses earlier than k(x). 



Substituting a = into fl27l) and fl28|) . we obtain 



h = r ^ l «M (48) 

/3(n-4)(n-3)(n-2)(n- l)^ 2 + 1 K ' 



and 



= - 5)( n - 4) 2 (n - 3) 2 (n - 2) 2 (ra - l)x 4 - 12(n - 4)(rt - 3)(rt - 2)x 2 + 3 

Similarly, formulas (129]) and (130]) lead to 



^ = ( n _ 2 y- 3 fc(vW 

dr [/3(n-4)(n-3)(n-2)(n- l)^ 2 + l] 3 



and 

fc(x) = i/( ra -8)(n-5)(n-4) 4 (n-3) 4 (n-2) 3 (n-l) ; 
9 I 

-2(n - 4) 3 (n - 3) 3 (n - 2) 2 (n - l) 2 (n + 13) (2n - 7)x 6 

+6(n - Af{n - 3) 2 (n - 2)(n - l)(n 2 - 8n + 45)x 4 

+18(n - 4)(n - 3)(2n 2 -n + 9)x 2 + 9} . (51) 

Fortunately, (3 dependence can be absorbed into the scaling factor in I and k. This is the 
reason why the analysis is relatively simple. 

To check the sign of h and dh/dr, we need to know the behavior of l(x) and k(x). These 
functions for n = 5 and n = 8 are plotted in FigJS] and FigJH respectively. In the case of 
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n=8 l(x) and k(x) 
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FIG. 4: The behavior of l(x) and k(x) in ten dimensions is shown. Apparently, l(x) has no positive 
root. Only k(x) crosses the x-axis. 



TABLE I: The lowest positive solution of l(x) = and k(x) = 0. 
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a 


b 


5 


jjg - 0.2041 


0.4076 


6 


0.1087 


0.1640 


7 

8 


0.09129 


0.09129 
0.05905 


9 




0.04171 



n = 8, as is shown in Figj4j l(x) is always positive. In fact, we can easily verify l(x) is always 
positive if n is larger than seven by calculating the discriminant of l(x) = 0. Otherwise, the 
equation l(x) = has a positive root. 

It is convenient to define a as the lowest positive solution of l(x) = and b as that of 
k(x) = 0. We see that l(x) is positive for < x < a and k(x) is positive for < x < b. 
The solutions a and b can be calculated numerically and the results are shown in table [H 
From this table, we see a is less than b for n = 5, 6 and 7. Note that the linear analysis 
is legitimate in seven, eight and nine dimensions if y/]3ipH < a,. Therefore, when the linear 
analysis is applicable, we have the relation \/J3i(jh < b. Thus, in the cases where the linear 
analysis is applicable, black holes turn out to be stable. In seven dimensions, for example, 
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comparing 

(USD and 1/2 V6, we can see that the linear analysis is good if \i is larger than 

32/5. 

For n > 8, is always positive and the solution (jSj) is stable if y//3ipH < b. This means 
that there exists a critical /i above which black holes are stable. 

To summarize this case, there is a critical \i in seven, eight and nine dimensions due to 
the limitation of the linear analysis and there also exists a critical \i in D > 10 due to the 
instability. 



D. a > and j3 > case in seven dimensions 



Now we will consider the third order Lovelock theory. Since the analysis is complicated, 
we discuss each dimension separately. From now on, we put (3 = to? (t > 0). From the 
perturbative point of view, it is natural to take the dimensionless parameter t to be at most 
order one. However, we keep it arbitrary parameter in our analysis. 

Now, formulas (J2ZD and ([28]) lead to 

24ta V + 12anf> + 1 v ; 

and 

l t {x) = 12(3 - 2t)x 2 + 6x + 1 . (53) 

As one can see, t dependence of l t (x) is a source of complication. Apparently, if t < 3/2, h 
is positive. Similarly, formulas (1211 and (I3"0l become 

dh = 2 W) 

dr (24ta 2 ip 2 + 12aip + 1) 3 1 } 

and 

kt{x) = 6912(3 - 2t)t 2 x 6 + 3456(4 - t)te 5 + 960t(6 + t)x A 

+ 2400tx 3 + 36(7 + 6t)x 2 + 42x + 1 . (55) 

From these expressions, we see the sign of l t (r) and k t (r) determines that of h and dh/dr, 
respectively. 

For < t < |, one can see lt(x) and kt(x) are always positive because we are considering 
the case aip > and then x > 0. This means that black holes are stable for any \i. For this 
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FIG. 5: The graphs of a{t) and b{t) are shown for t > 3/2 in D = 7. We clearly see the relation 
a(t) < b(t). 

natural choice of t, we do not have any instability under tensor perturbations. However, if 
we examine the scalar modes, it is likely that the instability shows up. 



lowest solution of lt(x) = and b(t) as that of k t (x) = 0. We numerically calculated a(t) 
and b(t) and the results are shown in Figj5l It is clear that lt{x) is positive for < x < a(t) 
and k t (x) is positive for < x < bit). As can be seen from Fig HI the relation a(t) < b(t) 
always holds. That means kt(x) is positive as long as the linear analysis is applicable, i.e., 
aipH < a{t). Hence, black holes are stable for these cases. Of course, there exist a critical 
mass for t > 3/2 due to the limitation of the linear analysis. 

E. a > and j3 > case in eight dimensions 

In seven dimensions, we have not seen the instability. However, in eight dimensions, we 
will see the instability of small black holes. 
From the formulas (I2T|) and (!28l) . we obtain 



For t > 4 



|, both k(x) = and kt(x) = have positive solutions. Let us define a(t) as the 




(56) 



120taV + 20aip + 1 ' 



where 



: t ( x ) = 960t 2 x 4 + 320te 3 + 24(5 - U)x 2 + I2x + 1 . 



(57) 
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15 16 17 18 19 20 21 22 23 24 25 
t 

FIG. 6: The behavior of a(t) and b(t) in D = 8 is plotted. Below i = 19.752, l t (x) = has no 
positive solution. 



From other formulas (j29]) and fl30l) . we have 

rfr (120toV 2 + 20a^ + 1) 3 ' K J 

where 

fc^z) = -2304000tV - 1536000tV + 16000t 2 (67 - 114t)x 6 + 24000t(13 - 10t)x 5 

+240(25 + 265t + 66t 2 )x 4 + 80(206t + 15)x 3 + 180(4 + 5t)x 2 + 66x + 1 . (59) 

Again, it is clear that the behavior of lt(x) and k t (x) determines the stability. 

First, we examine h{x). From the discriminant of Z t (x) = 0, we can see that lt{x) > in 
the region x > when t is less than 19.752. Next, we need to check the sign of kt(x). It 
is easy to show that k t (x) = has positive solutions for arbitrary t. We define b(t) as the 
smallest solution of k t (x) = 0. In FigJHl numerical results for a(t) and b(t) are plotted. From 
this figure, we can conclude that the solution (jHJ) for t < 19.752 is stable if aipH < b(t). In 
FigJTJ we numerically calculated the marginal /i, which satisfies aipn = b(t), as a function 
of t. A numerical fit gives us the formula for the marginal \i: 

/Wginai * (71.4671 t L47893 + 75.7542 t 2 - 61059 )a 5 / 2 . (60) 

Note that the error of this fitting is about one percent in the range t < 0.5. Therefore, black 
holes with fi less than this marginal value (160]) is unstable. It should be emphasized that the 
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FIG. 7: The marginal value of /x marg i na i divided by a 5 / 2 as a function of t is plotted with the solid 
line. A numerically fit to the data is plotted with the dashed line. 

larger the multipole orders i becomes, the shorter the time scale of the instability becomes. 
This is our main result in this paper. 

We find lt(x) = has two positive solutions when t is larger than 19.752. We define a(t) 
as the smaller one. The solution (jSJ) is stable for t > 19.752 if aipn < a(t). 

To conclude, for t < 19.752, we have found the instability of small black holes under 
tensor perturbations in eight dimensions. In any case, we have a critical mass below which 
black holes are unstable or the linear analysis can not be applicable. 

V. CONCLUSION 

We have studied the stability of static black holes in the third order Lovelock theory. 
We derived a master equation for tensor perturbations. Using the master equation, we 
have reproduced known results, which give a check of our master equation. The main 
purpose of this paper was the stability analysis of Lovelock black holes in seven and eight 
dimensions. We found there are cases where the linear analysis breaks down. If we restrict 
ourselves to the regime where the linear analysis is legitimate, black holes turns out to be 
always stable in seven dimensions. In particular, for a reasonable parameter t < 3/2, we 
found no critical mass. For t > 3/2, there is a critical mass due to the limitation of the 
linear analysis. In eight dimensions, for a reasonable range t < 1, we found a critical mass 
/^marginal ^ (71.4671 t 1 ' 47893 + 75.7542 t 2 ' 61059 )a 5 / 2 below which black holes are unstable. 
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Remarkably, the larger the multipole orders i becomes, the instability gets the stronger. It 
is interesting to note a similar result is found in cosmology in the presence of Gauss-Bonnet 
term 261 ] . 



We also examined a special case where the Gauss-Bonnet term disappears accidentally. 
In seven, eight, and nine dimensions, black holes are stable as long as the linear analysis 
is possible. However, if we regard the breakdown of the linear analysis as an indication of 
the instability, there always exists a critical mass. In more than nine dimensions, the linear 
analysis is always legitimate. We have found the instability for small black holes for these 
dimensions. Thus, in any case, we have a critical mass in dimensions 6 < D < 11. 

The results in seven and eight dimensions are similar to that of Einstein-Gauss-Bonnet 
theory. In Einstein-Gauss-Bonnet theory, there exists the tensor mode instability of small 
black holes in six dimensions and the scalar mode instability in five dimensions. Here, we 
should note that there exists no instability in eight dimensions in Einstein-Gauss-Bonnet 
theory. However, in eight dimensions, Einstein-Gauss-Bonnet theory is not a complete one 
in the sense of Lovelock. We found the tensor mode instability when we extend the theory 
to the third order Lovelock theory. It is not so absurd to imagine that, if we consider the 
scalar modes, there exists the instability in seven dimensions. If so, it is natural to expect 
this tendency continues. Thus, we conjecture that small black holes are unstable in any 
dimensions including ten dimensions. 

There are many remaining issues to be solved. First of all, we have to give the analysis 
of scalar and vector perturbations to complete our stability analysis. We also need to 
understand the meaning of the breakdown of the linear analysis. Moreover, it is worth 
investigating the fate of the instability. As the instability is stronger for higher multipole 
orders £, the resultant geometry would be weird. Finally, it is interesting to consider the 
Hawking radiation by taking into account the instability. The whole picture of the black 
holes at the LHC would be changed. 

Acknowledgments 

JS would like to thank Christos Charmousis for discussion on Lovelock black holes. JS 
is supported by the Japan-U.K. Research Cooperative Program, Grant-in-Aid for Scientific 



20 



Research Fund of the Ministry of Education, Science and Culture of Japan No. 18540262. 



[1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 
[arXiv:hep-ph/98033i~5] . 

[2] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 257 

(1998) [arXiv:hep-ph/9804398] . 
[3] C. Kokorelis, Nucl. Phys. B 677, 115 (2004) [arXiv:hep-th/0207234] . 

[4] S. B. Giddings and S. D. Thomas, Phys. Rev. D 65, 056010 (2002) [arXiv:hep-ph/0106219] ; 
S. B. Giddings and M. L. Mangano, Phys. Rev. D 78, 035009 (2008) [arXiv:0806.338Tl [hep- 
Ph]]- 

[5] H. Kodama and A. Ishibashi, Prog. Theor. Phys. 110, 701 (2003) |arXiv:hep-th/0305147| ; 
A. Ishibashi and H. Kodama, Prog. Theor. Phys. 110, 901 (2003) |arXiv:hep-th/0305185| ; 
H. Kodama and A. Ishibashi, Prog. Theor. Phys. Ill, 29 (2004) [arXiv:hep-th/0308128| . 
[6] R. A. Konoplya and A. Zhidenko, Nucl. Phys. B 777, 182 (2007) [arXiv:hep-th/0703231| . 
[7] K. Murata and J. Soda, Class. Quant. Grav. 25, 035006 (2008) |arXiv:0710.022T1 [hep-th]]. 
[8] M. Kimura, K. Murata, H. Ishihara and J. Soda, Phys. Rev. D 77, 064015 (2008) 

|arXiv:0712.4202l [hep-th]]. 
[9] H. Ishihara, M. Kimura, R. A. Konoplya, K. Murata, J. Soda and A. Zhidenko, Phys. Rev. D 
77, 084019 (2008) |arXiv:0802.0655l [hep-th]]. 
[10] K. Murata and J. Soda, Prog. Theor. Phys. 120, 561 (2008) [ arXiv:0803.1371l [hep-th]]. 
[11] B. M. N. Carter and I. P. Neupane, Phys. Rev. D 72, 043534 (2005) lj arXiv:gr-qc/0506103] . 
[12] H. K. Kunduri, J. Lucietti and H. S. Reall, Phys. Rev. D 74, 084021 (2006) 

[arXiv:hep-th/0606076] . 
[13] T. Oota and Y. Yasui. [arXw:0812.1623l [hep-th]. 

[14] H. Kodama, R. A. Konoplya and A. Zhidenko. larXiv:0812 .04451 [hep-th]. 
[15] D. G. Boulware and S. Deser, Phys. Rev. Lett. 55, 2656 (1985). 
[16] D. Lovelock, J. Math. Phys. 12 (1971) 498. 

[17] C. Charmousis, Lect. Notes Phys. 769, 299 (2009) [ arXiv:0805.0568l [gr-qc]]. 
[18] J. T. Wheeler, Nucl. Phys. B 268, 737 (1986); 
J. T. Wheeler, Nucl. Phys. B 273, 732 (1986). 



21 



[19] A. Barrau, J. Grain and S. O. Alexeyev, Phys. Lett. B 584, 114 (2004) [arXiv:hep-ph/0311238] . 
[20] G. Dotti and R. J. Gleiser, Class. Quant. Grav. 22, LI (2005) [arXiv:gr-qc/0409005] ; 

G. Dotti and R. J. Gleiser, Phys. Rev. D 72, 044018 (2005) |arXiv:gr-qc/0503117| . 
[21] I. P. Neupane, Phys. Rev. D 69, 084011 (2004) [arXiv:hep-th/0302132] . 
[22] R. J. Gleiser and G. Dotti, Phys. Rev. D 72, 124002 (2005) [arXiv:gr-qc/0510069] ; 

M. Beroiz, G. Dotti and R. J. Gleiser, Phys. Rev. D 76, 024012 (2007) |arXiv:hep-th/0703074| ; 

R. A. Konoplya and A. Zhidenko, Phys. Rev. D 77, 104004 (2008) |arXiv:0802.0267l [hep-th]]. 
[23] V. S. Rychkov, Phys. Rev. D 70, 044003 (2004) [arXiv:hep-ph/0401116] ; 

T. G. Rizzo, JHEP 0506, 079 (2005) [arXiv:hep-ph/0503163| ; 

T. G. Rizzo, Class. Quant. Grav. 23, 4263 (2006) |arXiv:hep-ph/0601029| . 
[24] R. C. Myers and J. Z. Simon, Phys. Rev. D 38, 2434 (1988). 

[25] M. Satoh, S. Kanno and J. Soda, Phys. Rev. D 77, 023526 (2008) |arXiv:0706.3585l [astro-ph]]. 
[26] S. Kawai, M. a. Sakagami and J. Soda, Phys. Lett. B 437, 284 (1998) |arXiv:gr-qc/9802033| . 



22 



